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Abstract 

We calculate the electronic polarizabihty in the superconducting state near 

extremum vectors Qq of the Fermi surface. A pole appears in the polarizabihty 

at frequencies oj near the superconducting gap 2A which leads to a sharp peak 

at LO just below 2 A in the lattice vibrations near Qq. A second order transition 

to a charge density wave state below the superconducting transition is shown 

to be unlikely. The results are compared with the recent inelastic neutron 

scattering measurements in compounds, ANi2B2C , yl = y or Lu. 
74.20.-z,74.25.kc 
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Strong momentum dependence of the electronic polarizability leads to characteristic fea- 
tures in the dispersion relation of phonons in simple metals-the Kohn anomalies. [|I| Tran- 
sition metals and compounds are particularly beset with wiggles and dips in their phonon 
dispersions. Strong electron-phonon interactions and local field effects are responsible to- 
gether with nesting features in the band structure for the anomalies and for the occurrence 
of charge density wave transitions 

It is well known that the leading order polarizability IIq{Q,uj), Fig. 1 (a), of an electron 
gas in d-dimensions for Q = 2kF for a spherical Fermi-surface or ellipsoidal Fermi-surface 
across the major and minor axes of the Fermi-surface in the normal state is 

ReU^{2kF,uj) ~ A^(0)[1 + {N{0)ujy~Hn{—)]. (1) 

We show the existence of singularities in the polarizability at ^ 2A in s-wave super- 
conductors near the extremum vectors of the Fermi-surface which are stronger than those of 
the normal state at ^ 0. By an extremum vector we mean Qq such that ^"^^ 
are both on the Fermi-surface and for small deviations \k\ << \Qq\ 

P 

Here k\\ is the component of k parallel to Qo and ki± are the components perpendicular to 
Qo where i runs over d — 1 values. For example, Qq are the major and minor axes for an 
ellipsoidal Fermi- surface. Qq is also a nesting vector when mi± for one or more i is infinity. 

These singularities strongly affect the phonon spectrum of superconductors under certain 
conditions. We explain the remarkable difference between the phonon spectrum of YNi2B2C 
in the superconducting state and that of the normal state recently observed by Kawano et 
al 0] and in LuNi2B2C by Stassis et al 0. In the normal phase of {Y, Lu)Ni2B2C the 
transverse acoustic and optic dispersion curve in the TZ direction has a strong temperature 
dependent dip around Q = (0.5,0,8). In the Y compound, below the superconducting 
transition temperature, the spectral weight splits into two parts, a very sharp resolution 



limited line at an energy ~ 2A and a broad part around the plionon peak of the normal 
state with the total spectral weight conserved. Weight is transferred from the sharp peak to 
the broad peak at low temperatures by a magnetic field. 

As a model, we consider a simple ellipsoidal Fermi-surface, = Y^i and calculate 
Ii{Q,uj) in the superconducting state for Q near the extremum vectors Qq. 

We write the model Hamiltonian which describes a system of electron interacting via a 
potential V, 

in Nambu notation 0, where the electron and annihilation operators are written as two- 
component vectors, 



and the r's are Pauli matrices. The potential V{k, k', Q) includes the electron-electron 
Coulomb repulsion as well as the attractive interaction mediated by the phonons. As usual 
we rewrite Eq. (3) as 

H,i = Ho + H, (5) 
where Hq is the BCS reduced Hamiltonian, 

^o = E^l(efc^3 + ArOvl>,^, (6) 

k 

and Hi now includes the second term of Eq. (3) minus the electron-electron effective inter- 
action absorbed in Eq. (6). 

The lowest order polarizability. Fig. 1 (a), at T = is 

no(Q,^) = -ij -j^Tr[T,G{k + Q,uj + uj')T,G{k,uj')]. (7) 

where G is the single particle Green's function for the BCS reduced Hamiltonian. 



k 



where Ei = (e^ + A^) 



To evaluate Eq. (7), it is helpful, after performing the frequency integral, to divide 
the momentum space into parts in which and -Eg^g are close to A, say between A and 
sA, s ~ 0(2), and those in which they are larger than sA. In the latter E^j:_^0 may be 
replaced by eg^+g, so that the integrand is the same as in the normal state. Any new 
features expected in the superconducting state can arise only from the singular density of 
one-particle states for energies near A, i.e., from regions of momentum space where both Ej: 
and E^^^ are near A. The contribution of this region in three dimensions is, 

Reng^(Q,cu) ~ ^^^^^sgn{6uj) In 



r 6u <o 

where 6uj = lo — 2A(1 + A^(0) = " ^ where v is the volume of the imit 

cell, and q defined to be {Q — Qo)\\ satisfies qC, « 1. 

In two dimensions, IleIlQ'^{Qo,uj) is proportional to |2A — uj\^^ while in one dimension 
(where the treatment here is quite inadequate) it is proportional to |2A — lu\^2. For non- 
extremum Q spanning the Fermi-surface, IieIlQ'^{Q,uj) is proportional to ln|2A — uj\ in two 
dimensions and |2A — co]^ in three dimensions. 

The contribution of the other regions differs from the normal state only through a lower 
cut off at sA so that \n{EF/uj) in Eq. (1) becomes \n[EF/{uJ + sA)]. Since we will be 
concerned with u ^ 0(2 A), this change is not of consequence, so that 

Reno(Q, uj) - ReU'o'iQ, uj) + ReW^iQ, uj) (10) 

In the random phase approximation. Fig. 1 (b), 

n(g, u) = UoiQ, + V{Q)Uo{Q, uj)), (11) 
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The direct electron-electron interaction is more important in V{Q) in Eq. (11) compared to 
the phonon- induced interaction because of the cut-off Ep in the former compared to Debye 
frequency in the latter. Vertex correction to Eq. (11) only change the magnitudes calculated 
here without affecting the singularities. This is unlike the case Q — > 0. P| 

Since IIq'^{Q,lj) is singular for Q = Qo and u = 2A while IIq{Q,u) is smooth and much 
larger than IIq'^ for all u) except near 2A, polarizability is approximately given by 

Ren(Qo, ^ Ren'"(Qo, uj) u ^ 2A 

Ren(Qo,w) ^ Ren"(go,c^) elsewhere, (12) 

Here Il^{Qo,u!) is the normal state polarizability and 

Reff'=(Qo, ^) ^ zl^i^ . (13) 

(3 + VReUf,%Qo,uj) 

In Eq. (11), we have used RcIIq ^ — A^(0) and /3 ~ 1 — VN{0). In three dimensions W'^{q, to) 
always has a pole at 

a; = 2A(l-e"^), (14) 
where 7 = AN{0)V/ Ep « 1. Note that this pole is very close to 2A and has weight 

r ~ 2EFe'^/V « 1. (15) 
The phonons couple to the electronic charge density, 

Hei-,H = E 9Q{bQ + bl^)^^.rs^^„ (16) 

k 

where gq is the electron-phonon coupling constant. The phonon spectrum in a metal is 
described by 

D-\Q, uj) = D,\Q, u) - S(g, u), (17) 

where T,{Q,uj) = gQll{Q,uj) and Do{Q,uj) contains no renormalization due to electron- 
phonon interactions. Using Eqs. (12) and (13) in Eq. (17), 
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D{Qo, uj) = [u- UQ, - + zTQ.eiu - 2A)]-' 

2 

where and Tq^ are the phonon frequencies and hnewidths of the normal state obtained 
from the normal contribution IP{Qo,uj) to II{Qq,uj). The 6'-function in Eq. (18) takes into 
account that Imn(Qo;^) = for u < 2A. Actually, we should include in Eq. (18) also a 
contribution to Tqq due to ImW^^Qo, u) which increases damping of the normal state if cuqo 
is between 2 A and ^ 0(4A). 

The electron-phonon coupling function Qq^^ ~ XCoEp where ui is the order of the Debye 
frequency and A is the dimensionless electron-phonon coupling constant, typically of 0(1). 
Then {gq^rY^'^ is typically much less than 2A. In that case and for ujq^^ > 2A, Eq. (18) 
predicts a pole at 

^ ^ 2A - Igl,r/{UQ, - 2A) ^ z., (19) 
and peak with width Fq^ at 

1 



UJ 



2 

The spectral function of the pole near 2A is 



^Qo + o^?Q.//(^Qo-2A). (20) 
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with the rest in the peak near uq^. If uq^ is close to 2 A, the pole is at 2 A — gq^r^^"^ and 
shares weight equally with the peak. The spectral function, S{Qq,uj) = —^lmD{Qo,u) is 
illustrated in Fig. 2. For g 7^ 0, the new mode increases in energy by 0(A^^g^/4), and is 
over damped for q > 

In comparing with experimental results in {Lu, Y)Ni2B2C, one should bear in mind 



that the d = 3 results in normal state are explained in a detailed calculation |[Tl| of HoiQ, ^) 
based on the actual complicated band structure and phonon renormalizations using electron- 
phonon coupling incorporating local filed effects given by methods of Ref. p|. We have 



studied only the effects due to superconductivity assuming ojq^^ in Eq. (18) to be the observed 
normal state value just above T^. The band structure of {Y, Lu)Ni2B2C shows the extremum 
vectors in the (a*, h*) plane for Qq ~ (0.5, 0, 0), and weak dependence on momentum in the 
c* direction. The local field effects quite generally enhance the magnitude of the effects 
in the superconducting state just as they do in the normal state. But since they are not 
included in our calculation not too much reliance can be put on the estimate of magnitude 
of the coefficients in Eqs. (19)-(21). Moreover band structure calculations shows significant 
nesting near Qq of the anomaly. For well nested surfaces, the d = 2 result quoted below Eq. 
(9) is relevant and the effects we discuss are correspondingly stronger. 

YNi2B2C is a particularly good candidate for the effects discussed here because the 
transverse acoustic phonon at Qo — (0.525, 0, 8) softens to an energy ^ 7 meV just above 
Tc, not too far above 2A ~ 4.3 meV. Thus the weight transferred to the peak near 2A is 
significant and estimated from Eq. (21) to be the experimental value 0.3 if r/2A = 

0.4A^(0) and 7 ~ 0(0.01) at T = 0. Actually, we have no independent way of estimating r. 
This value is reasonable only if /? ^ 0(0. 1), which implies a phonon close to an instability 
in the normal phase, as in these compounds. 

A simple extension of the calculation above shows that as the temperature is decreased 
and A(T) grows, weight is transferred to the sharp peak at the expense of the broad peak 
with the total weight constant. Similarly a magnetic field depresses the weight and the 
frequency of the sharp peak. These are as observed [Q. The theory would also predict a 
temperature dependence to the frequency of the peak ^ 2A(T). In the experiment 
the frequency of the new peak appears not to follow the BCS dependence of 2A(T). This 
may be so because of strong-coupling effects. Note also that the relative intensity of the 
peak drops rapidly as its frequency decreases below 2A(0), i.e., near T^, making it harder to 
observe. 

We finally comment on other possible explanation of the data. The observed peaks 
cannot be the large Q-extensions of the amplitude modes of superconductors near Q — 
observed in Raman scattering [|l^ since such modes are heavily damped for \Q\> 
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Another possibility is that the peak is just the narrowing of the hnewidth as the phonon 
frequency goes below 2A. This would predict that the number of peaks for a fixed Qo as 
a function of u remains the same in the normal state and the superconducting state. In 
Kawano et al's data 0], an extra new peak for T < Tc clearly appears, as in the theory here. 
In the data of Stassis et al in LuNi2B2C 0, an extra peak is not visible. In this compound 
the phonon frequency of the soft transverse acoustic mode in the normal state just above Tc 
is lower than that in y-based compound, and in fact close to the sharp peak appearing well 
below Tc. Recall that if ojq^ is close to 2A(0), the two peaks in the superconducting state 
are split only by gQ^r^^"^ as discussed below Eq. (21). The two peaks may then lie within 
the experimental resolution of about 1 meV. Further experiments are required to elucidate 
the situation. 

An interesting conclusion from our work is that a second order transition to a charge 
density wave state at a temperature, Tcdw, below the superconducting transition tempera- 
ture, Tc, may not be allowed. A second order transition requires UQ^iT) ^ as T ^ Tqdw- 
But co'Q(,(T) is bounded from below by 2A(T) for Tc > Tqdw, because of the singular con- 
tribution to the phonon self-energy as cuq^ approaches 2A(T). The variation of the phonon 
frequencies as a function of temperature is schematically illustrated in Fig. 3. A first order 
transition is, of course, allowed. These matters need further investigation. 

In summary, we find the singularities in the polarizability of s-wave superconductor at 
nesting vectors in three and two dimensions, which manifest themselves as poles in the 
lattice response function at u just below 2 A with its weight dependent on the frequency 
of the normal state phonon at the nesting vector. These results are in general accord 
with observations. A straightforward extension of these results can be made to d-wave 
superconductors. If the maximum of the gap function in such superconductors is along the 
nesting direction, the results for d-wave superconductors are expected to be similar to those 
obtained here. 

We are grateful to Dr. H. Kawano and Prof. C. Stassis for discussions of the experimental 
results and to Prof. B. Altshuler for pointing out an error in the calculation. We have also 
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learned that Dr. S. K. Sinha(private communication) has arrived at conclusions similar to 
those in the present paper. One of authors (H. Y. K.) thanks E. Abrahams, S.-W. Cheong, 
Y. B. Kim, and P. B. Littlewood for helpful discussions. This work is supported in part by 
the Korea Research Foundation and NSF Grant No. DMR-96-32294 (H. Y. K.). 
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Figure Captions 

Fig. 1: (a) Lowest order of polarizability IIq{Q,uj), and (b) Polarizability Il{Q,u!) in random 
phase approximation. 

Fig. 2: Phonon spectral function S{Qo,uj) calculated from Eq. (21) for various uq^ where 
<^Qo = ^Qo/2A, and r/2A = 0.4Ar(0). 

Fig. 3: Schematic behavior of the phonons as a function of temperature. For Tc below the 
charge density wave transition temperature Tcdw, there is only a single branch luq^ — > as 
T — > TcDW: as shown by the dashed line. For Tc > Tcdw, the soft phonon branch sphts into 
two for T <Tc, the lower one uji is slightly below 2A(T) and the upper one, u;2 considerably 
broadened with its lower edge at slightly above 2A(T). 
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Fig. 1 (b) 
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